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Abstract 

Let K he a random variable following a truncated exponential dis¬ 
tribution. Such distributions are described by a single parameter here 
denoted by 7 . The determination of ^ by Maximum Likelihood methods 
leads to a transcendental equation. We note that this can he solved in 
terms of the inverse Langevin function. We develop approximations to 
this guided by work of Suehrcke and McCormick. 
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1 Introduction, applications 

Preamble, not for journal publication. This paper was in 2007 submitted to 
Solar Energy as it was a slight improvement on the theory part of a paper 
in that journal. As it lacked observational data, that journal suggested it 
wasn’t appropriate as it stood. (WRF simulations have been considered, 
but not implemented.) In 2012 another journal proposed a Special Issue 
on Special Functions and Probablity and Statistics Applications, but this 
proposal didn’t eventuate. The arXiv version here is the 2012 version of the 
paper. 


The frequency distribution of insolation values is important for the pre¬ 
diction of the performance of systems involving solar heating. The fractional 
time distribution for the clearness A, for K rr^\r^ < K < Amax can be modelled 
by 


F{K) 


exp(7Amin) - exp(7A) 
exp(7Amin) - exp( 7 Amax) ’ 


( 1 ) 
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see [2] . (Of course, F{K) = 0 for K < -fCmin and F{K) = 1 for > Kmax-) 
This is the doubly-truncated exponential distribution. We remark that 7 
could be negative. Equation ([T]) appears in equation (3.7.4) of [9] p72 and 
equation ( 1 ) of [ 11 ]. ETmax and iCmin are the maximum and minimum of the 
daily clearness index over the period of interest (and a common approxima¬ 
tion is that = 0.05). Simulations of solar water and air heaters are 

discussed in |9]. Oher situations - and one which the author encounted is 
described in |5] - might also be simulated in like manner. 

The truncated exponential distribution has many other applications in¬ 
cluding distributions of earthquakes, of forest-fire sizes, raindrop sizes, re¬ 
liability modelling, etc. See mmm for example (noting that, in some 
applications, Tfmax and 7 are all unknowns). The Langevin function, 

defined in ^ is easily recognized in formulae in these papers. For example, 
equation (2.7), h 2 {i') = r, in [3], is L{v/2) = 2r — 1. 

If one supposes ETmax and iCmin values to be known, it remains to esti¬ 
mate 7 . To do this, first calculate 



The quantity 7 can be estimated by solving K = fix { 7 ), where K is the 
mean of the observations. See [9] p72 equation (3.7.5) and [H] equation 
(2). The 7 so found is the method of moments estimate and, in fact, the 
maximum likelihood estimate of 7 . 

To shorten the formulae which follow, define a (an average) and 6 (half 
the difference) by 


— 2 ('^™Eix T T^min), <5 — 2^^' 


max 


Using Eimin = a —6, and it'max = a-|-5 in the right-hand side of equation ([2|) 


= - \- a + 6 coth(75) 


(3) 


so equation ([ 2 ^) becomes 



(4) 
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2 The Langevin function 


The purpose of this genuinely elementary note is to call attention to the 
fact that the function occuring on the right-hand side of equation Q is the 
Langevin function, and various consequences of this. For example, 7 can 
be expressed exactly in terms of the inverse Langevin function. Both the 
Langevin function and its inverse have been studied widely in other contexts: 
see, for example, HEmn]. The Langevin function is defined, with x 7 ^ 0, 
by 

, , , , 1 d , /sinh(x)\ ,,, 

L{x) := coth(x)-= — log ( - ) (5) 

X ax \ X J 

and, for x = 0, by L(0) = 0. The function L is a monotonic increasing 
function mapping the real line into the interval from —1 to 1 , and it satisfies 

= 1 _ (6) 

ax X 

The inverse Langevin function is the function L~^ inverse to L. Of course, 
the computation of L~^ remains a numerical issue, just as the solution for 
7 of equation m) is. However, with the knowledge of the name of the 
function, the solar applications can share algorithms and code with the 
other applications, if it is so desired. Knowing the name of the function 
involved assists in the search for numerical software. Numerical codes, in C, 
are available for the inverse Langevin function. 

With 


X — .^min)) V — 


{K - a) 2K - (Kmax + Kmm) 


K — K 


equation dH) can be written y = L(x). Consequently 


7 = \y)- 


(7) 


Suehrcke et al. [m gave an approximation to 7 solving equation (j4]) 
which amounts to approximating x = L~^{y) by 

X « ^tan(-^). 

Suehrcke used a polynomial expression for A, finding its coefficients in the 
best way to be useful over the range of values of iLminj TLmax and K occuring 
in the clearness application. Restricting to the one tangent evaluation, this 
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can be used in other ways to approximate L~^. Let r = tan( 7 r?// 2 ) = 
tan( 7 rL(x)/ 2 ). By fitting for small r and for large r using Maple we find 
that 


X 


IT 


-y 


/1 + bnr'^ \ 


where 



20tt^ - 144 
7r2(60 — Stt^) 


0.508. 


This approximates L~^ with a relative error of about 0.3% over the entire 
range and even smaller relative errors can be obtained by using higher degree 
polynomials in the numerator and denominator. 

The Langevin function - and nothing worse than this - arises in the 
formulae for all the higher moments and for all the higher cumulants in the 
truncated exponential distribution. The moment generating function M, 
and the cumulant generation function C, are defined by 


M{s) := [ 
Jh 

One easily finds 


exp(siL)—— dK, 


dK 


C{s) := log(M(s)). 


C{s) = sa + 


, / sinhfx + s5) \ 


, / sinh(x) \ 


Differentiating this gives 

dC 

C'{s) = — = a + 6L{x + S(5). (8) 

ds 

The fact that C"(0) = Hk{7) recovers equation ([3|). Equations ([ 6 ]) and ([ 8 |) 
show that each higher derivative of C can be expressed as a polynomial in 
L = L{x + s6}, S, and l/(x + sS). The most important cumulant after the 
mean fiRij) is the second, the variance yi'y)- 

9 T 

alh) = C"{0)=5\l-y-—), (9) 

where L = L{'y5) = — a)/5. (The right hand side can be expressed, 

equally well, in terms of coth(x) as in terms of L{x).) 

There is a practical use for equation Q. Let the first and second cu¬ 
mulants from the data be denoted K (for the sample mean) and s\ (for 
the sample variance). Next, denote by v the value of the right-hand side 
of equation ([9]), when 7 is given by that solving K = and 
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is replaced by K. would give some indication of how well the truncated 
exponential distribution (with the given Kmin, -fi^max) fits the data. 

There is a more general context for some of this, for example, in the 
context of natural exponential families with quadratic variance functions. 
Equation ([9]) is an example of a ‘quadratic variance function’. The untrun¬ 
cated versions of these distributions are treated in [7] , but further considera¬ 
tion of such matters takes us too far from the Langevin function observation 
which is the subject of this note. 

3 Conclusions 

We hope that practitioners in the different disciplines using truncated ex¬ 
ponential distributions will be able to benefit from sharing experiences. In 
particular, in this paper, we have called attention to the Langevin func¬ 
tion, and the existing methods of approximating its inverse appropriate to 
estimating 7 . 
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